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*~ 5 Abstract 

Let N(t) denote the eigenvalue counting function of the Lapla- 
cian on a compact surface of constant nonnegative curvature, with or 
without boundary. We define a refined asymptotic formula N(t) = 
At + Bt l l 2 + C, where the constants are expressed in terms of the ge- 
ometry of the surface and its boundary, and consider the average error 
A(t) = \ J * D(s) ds for D(t) = N(t) - N(t). We present a conjecture 
for the asymptotic behavior of A(t), and study some examples that 
support the conjecture. 
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>• "The mills of God grind slowly, yet they grind exceeding small." 
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1 Introduction 

m 

For any positive self- adjoint operator with discrete spectrum 
X < Ai < A 2 < A 3 < • • • — > oo (repeated according to multiplicity) 

we consider the eigenvalue counting function 

N(t) = #{A, < t}. (1.1) 
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Often there is a predicted asymptotic approximation N(t), and the problem 
arises to estimate the error 



D(t) = N(t) - N(t). (1.2) 
Here we investigate some examples where the average error 

1 f* 

A(t) = - / D(s)ds (1.3) 
t Jo 

is better behaved, following up on [JS]. In our examples we deal with — A on 
a compact surface S which is either flat or has constant positive curvature. 
The Weyl asymptotic formula has N(t) = At + o(t) where the constant A 
is given by A = ^Area(S'). If S has a reasonable boundary then there is 
a second term in the asymptotics, N(t) = At + Bt 1 ^ 2 + o(t 1 / 2 ), where the 
constant B depends on the boundary conditions. If we impose Neumann 
boundary conditions on all of dS, then B = j-length(dS'), and if we impose 
Dirichlet boundary conditions then B = — ^-\ength(dS). We may also con- 
sider mixed boundary conditions, splitting dS = 8Sn U 8Sd and imposing 
Neumann boundary conditions on dS^ and Dirichlet boundary conditions 
on dSr>. In that case B = ^length((9S'Ar) — ^length(9S , £i). 

For our purposes we will need a more refined asymptotics with an addi- 
tional constant term: 

N(t) = At + Bt 1/2 + C (1.4) 

where A and B above and C is given as follows. We assume that dS 

is piecewise smooth, with a finite number of corners with angles {9j}- We 
will write C = C\ + C 2 + C 3 , where C\ is the contribution from the corners, 
C 2 is the contribution from the curvature of the smooth arcs in dS, and C 3 
is the contribution from the curvature of S. Thus C 3 will be zero if S is flat, 
and 

Cs = ^K 2 (S), (1.5) 

where K 2 (S) is the total curvature, the integral of the curvature over S, in 
the general case. Thus for any sphere K(S) = 4n and C 3 = |. For C 2 we 
will take the integral of the curvature K\ of the boundary, 



C 2 = Azf K i- (1-6) 
las 
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Here the curvature is taken with respect to S, so it will be multiplied by — 1 
on the interior portions of the boundary. If dS has no corners and consists 
of an outer boundary and N inner boundary curves, then C2 = 

To describe the constant C\ we define 




If we impose Neumann or Dirichlet conditions throughout, then we take 

where the sum is over all corner points: For mixed boundary condition we 
subdivide the boundary at corner points, and impose Neumann conditions 
on the arcs between some corner points, and Dirichlet conditions on the 
remaining arcs. We then sort the corner points into those {9'j} where the 
same boundary condition is imposed on both incident arcs, and {9"} where 
different boundary conditions are imposed. Then we take 

ci = + E W) - W))- (i-9) 

We may also allow the surface to have a finite number of cone point 
singularities with cone angles {%}, < a j < 2n. In that case we add to C\ 
the value 

E 2 ^(f)- ( L1 °) 

Conjecture, (a) Suppose S is flat. Then there exists a uniformly almost 
periodic function g of mean value zero, such that 

A(t) = (?(t 1/2 )r 1/4 + 0(r 1/2 ), (1.11) 

and more generally there exists a sequence {gj} of uniformly almost periodic 
functions, with g\ = g, such that for all N 

N 

A(t) = J29j(t 1/2 )t~ j/4 + 0{r^/ 4 ) (1.12) 

3=1 

(b) Suppose S has constant positive curvature. Then there exists a uniformly 
almost periodic function g of mean value zero, such that 

Ait) = g(yjt + \) + Or 1 / 2 ), (1.13) 

3 



and more generally there exists a sequence {gj} of uniformly almost periodic 
functions, with g\ = g, such that for all N 

m = E*(\A+I )r( " 1)/2 + 0{rN/2) (L14) 

3=1 

In this paper we present a number of examples that support the conjec- 
ture. Of course this is exactly backwards, since the examples were worked 
out first and the conjecture was concocted to agree with the examples. All of 
the examples are sufaces for which it is possible to compute the spectrum of 
the Laplacian exactly in terms of trigonometric polynomials for flat surfaces 
and spherical harmonics for positively curved surfaces. In particular, they 
are highly symmetric. So they provide only weak evidence for the conjec- 
ture, and the methods of this paper do not provide a pathway to attack the 
conjecture. Nevertheless, the conjecture seems interesting enough that it is 
worth submitting to the mathematical community to stimulate further work. 
We should also point out that although all the almost periodic functions in 
our examples of positively curved surfaces are in fact periodic, we have not 
made periodicity part of conjecture (b). 

For surefaces without boundary the following relationship appears to 
hold: the set of frequencies in the trigonometric expansions of the almost 
periodic functions coincides with the set of lengths of closed geodesies on the 
surface. At present we have no explanation for why this should be the case. 

There should be an analogous conjecture for surfaces of constant nega- 
tive curvature. The expression for the refined asymptotics N(t) should be 
the same, but it is not clear what should replace (1.11) and (1.13). There is 
a vast literature on the spectrum of the Laplacian on compact Riemann sur- 
faces, or even surfaces of finite volume, and the relationship with the lengths 
of closed geodesies, going back to the Selberg trace formula (see [Bu], [Sa] 
for example). Although the spectral projection operator on the hyperbolic 
plane is known explicitly ([Ta], [S]), there is apparently no known exact com- 
putation on any quotient space. Since we have no examples to compute, we 
will not venture a conjecture. 

There is also a vast literature on the asymptotics of the trace of the heat 
kernel ([G]). The trace of the heat kernel is given as a smoothing of the 
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eigenvalue counting function 



h(t) = J^ 6 ^ = 1 \ N(s)e- St ds, 



(1.15) 



so asymptotics of N(t) translate immediately into asymptotics of the heat 
kernel trace. If we substitute N(s)—N(s) = -^(sA(s)) in (1.15) and integrate 
by parts we obtain 



so asymptotics for A(t) also translate into asymptotics of h(t). To go in 
the reverse direction requires the application of a Tauberian Theorem, and 
cannot reveal the refined asymptotic statements that we are interested in. 
However, we can check that the refined asymptotics N(t) that we use are 
consistent with the known asymptotics for h(t) ([BS], [K], [G]). Our average 
function A(t) is a much cruder smoothing of N(t) than the trace of the heat 
kernel. This is a double edged sword. On the one hand, we can't expect as 
nice behavior. On the other hand, the rougher behavior reveals some inter- 
esting new features; for example, the almost periodic functions. 

If the boundary is smooth, we observed that the contribution C2 to the 
constant term is given by "topological data," namely C 2 = ^ where g is the 
genus of the surface. The analogous statement for the trace of the heat kernel 
is the grand finale to the famous paper of Mark Kac [K]. It is interesting to 
note that Kac's argument is based on approximating the smooth boundary 
by a polynomial curve, and he is able to say something about the polyno- 
mial case in terms of some frightful integrals, but only in the limit does the 
result become comprehensible. (An explicit formula for the polynomial case 
is given in [BS].) 

In our approach, the polygonal case is quite explicit, but we may also 
check that there is continuity in approximating the smooth boundary by a 
polygonal boundary. To simplify the discussion, assume that the surface is 
a simply connected convex planar domain and has a smooth outer boundary 
dS. Subdivide the boundary at n points {xj}, and connect them by line 




where 



(1.16) 




(1.17) 
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segments to form a polygon P. Then the constant term in N for the polygonal 
domain is entirely given by C\, namely (1.8) where 9j is the angle at Xj. As 
n gets large the values of Oj approach n from below, and 

24^(0,-) = 2(1 - % + 0((1 - ^) 2 ). (1.18) 

7T 7T 

A simple geometric argument shows 

J20j = (n-2)n (1.19) 

and so 

Ci = £mo = ^( 2 ^ - £ E **) + * = I + * ( L2 °) 

where the remainder i?isO(^(l — ^-) 2 ) and tends to zero in the limit. Thus 
the constant term for the polygonal approximations approximates the con- 
stant term (C 2 ) for the smooth domain. It is straightforward to localize this 
argument to polygonal approximations to portions of the boundary of any 
surface. 

We want to emphasize that it is important to average the error D{t) in 
order to get the asymptotic behavior. In all our examples the function D(t) 
is unbounded. It is the fact that it is positive and negative that allows can- 
cellation to produce bounded behavior for A(t). Compared with the values of 
D(t), the constant term in the asymptotics N(t) is indeed "exceeding small." 
The fact that it nevertheless shows up in the average seems truly remarkable. 

Similar results should be valid in higher dimensions, but they will require 
a different type of average. We leave this to the future. 

We now outline the examples that take up the rest of this paper. There 
are two basic examples, the torus discussed in [JS] and the sphere, discussed 
in section 4. In all the other examples, the eigenfunctions on the surface may 
be extended to eigenfunctions on a torus or a sphere, so the spectrum on 
the surface is a subset of the spectrum for one of our basic examples. Our 
task is then to identify exactly the subset, and show how the conjecture for 
the basic example yields the conjecture for the surface. This requires only 
elementary reasoning, but the arguments are a bit technical. Of course we 
have to be very careful, since the value of the constant depends on getting 



6 



exact statements. It may seem that we are working out a lot of examples 
using very similar arguments. However, we found that we needed all the 
examples to help formulate and confirm the conjecture. We have tried to 
present enough detail so that the reader can verify the correctness of the re- 
sult, without excessively repeating the framework of the reasoning. We have 
used what we hope is self explanatory notation, but it changes from example 
to example. So, for example, we write N(t) for the counting function for 
the surface under discussion (or Niy(t) or No(t) if there is a boundary with 
Neumann or Dirichlet boundary conditions). If we have to recall a counting 
function from a previous example we will write N T (t) for a torus T, etc. 

In section 2 we discuss examples that are polygons in the plane with either 
Neumann or Dirichlet boundary conditions throughout. These examples are 
arbitrary rectangles, and certain special triangles: right isosceles, equilateral, 
and 30° - 60° - 90°. 

In section 3 we return to the same surfaces, but deal with mixed bound- 
ary conditions, Neumann on part of the boundary and Dirichlet on the re- 
mainder. We are able to handle all possibilities for the rectangle and right 
isosceles triangle, but none for the equilateral triangle and only some for the 
30° — 60° — 90° triangle. In this section we also consider an arbitrary cylinder 
and Mobius band. 

In section 4 we consider the sphere, the hemisphere and the projective 
sphere. In section 5 we consider lunes and half-lunes. These examples are 
very useful because they provide a wider variety of corner angles than the 
previous examples. 

In section 6 we discuss surfaces with point singularities. These include 
the flat projective plane discussed in [JS], the surface of a regular tetrahedron 
discussed in [GKS] , half tetrahedra, and glued-lunes. 

In section 7 we re-examine some of our examples of surfaces that have 
a finite group G of isometries. The question is how the spectrum sorts into 
the eigenfunctions with prescribed symmetry, given by the irreducible rep- 
resentations {-kj} of G. A heuristic suggested in [S] is that the proportions 
Nj(t)/N(t), where N t (t) is the counting function for eigenfunctions of ttj 
symmetry types is asymptotically (dimiTj) 2 /#G. Here we work out refined 
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asymptotics Nj(t) for Nj(t), by reducing the computation of Nj(t) to previ- 
ously worked out examples for a fundamental domain of the G action. The 
leading term At is as predicted, but the next term Bt 1 ^ 2 may be positive, 
negative, or zero, showing that the individual representations are somewhat 
overrepresented or underepresented in the spectrum. Because we only have 
a few examples, we are not able to formulate a conjecture for the behavior 
in the general case. This is another interesting open problem for the future. 



2 Polygons 

Consider the rectangle R of side length, a, b, and let T denote the torus 
of side length 2a, 26. Any Neumann eigenfunction on R extends by even 
reflection to an eigenfunction on T, and similarly a Dirichlet eigenfunction 
extends by odd reflection. In either case we obtain roughly a quarter of the 
eigenfunction on T, but in fact we can give a precise formula relating the 
counting functions N N and N D for R with N T for T. 

The eigenfunctions on T have the form 

e(j,k) = e* i & + *v\ j,keZ, (2.1) 
with eigenvalue + E ^f-- The Neumann eigenfunctions on R have the form 

3 ^ 

c(j,k) = cos t\-x cos TT-y, j > 0, k > 0, (2.2) 
a b 

while the Dirichlet eigenfunctions have the form 

7 A/ 

s(j,k) = sin Ti-x sin 7r— 2/, j > 0, > 0, (2.3) 
with the same eigenvalue. 



Lemma 2.1. W^e /iawe 



N N (t) = -N T (t) + 



iV D (0 = - A N T (t) 



1 


"at 1 / 2 " 


1 


'bt 1 ' 2 ' 


2 


7T 


+ 2 


71 


1 




1 


rot 1 / 2 ! 


2 


7T 


~~ 2 


7T 



+ 



1 

4 



(2.4) 
(2.5) 
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Proof. In the generic case j > and k > 0, we have four eigenfunctions 
e(±j, ±k) contributing to Nt for one eigenfunction c(j, k) or s(j, k) con- 
tributing to Nn or A^d, giving rise to the \Nxit) terms in (2.4) and (2.5). 
We then have to correct for the nongeneric cases. Note that e(0, 0) con- 
tributes once to N T and N N but not to N D , so this gives rise to the constant 

2 - 2 

terms. When k — and j > 0, we have e(±j, 0) for z y~ < t contributing to 
N T (t) and c(j, 0) contributing to N N (t) but not N D (t), so this gives rise to 

the term ±| in (2.4) and (2.5). Similarly the case k > and j = 

gives rise to the term ±£ . □ 



Note that the function [x] is on average a; — |, so we may rewrite (2.4) 
and (2.5) as 



N N (t) = ±N T (t) + \ 



N D (t) = -N T (t) 



at 1 ' 2 



7Y 



at 1 ' 2 



TV 



+ 



+ 



+ 



bt 1 ' 2 



7Y 



bt 1 ' 2 



7Y 



1 

4 



+ 



+ 



(2.6) 



(2.7) 



Now the constant term is the same in both equations. We define the refined 
asymptotics 

2a + 2b 



N„(t) = -t + 



N D (t) 



ab 



-t 



4ir 
2a + 26 



t^ 2 + 



f l/2 



(2.8) 
(2.9) 



A-K 47T 

Note that ab is the area of R and 2a + 2b is the length of the perimeter of R 
as predicted, and \ = ^ + ^ + ~[Q J r-[Q as > predicted. As usual we define the 
discrepancy D(t) and average A(t) in each of the three cases. 



Theorem 2.2. Both A N and A D satisfy 

A(t) = g(t^ 2 )t-^ 4 + Oit- 1 ' 2 ) ast^oo 
where g is an almost periodic function of mean value zero. 
Proof. We have 



(2.10) 



D N (t) 
D D (t) 



±- 



at 1 ' 2 



+ 



at 1 / 2 



TT 



±- 



bt 1 ' 2 



TT 



1 

2 



bt 1 ' 2 



TT 



1 



D T (t). 



4 

(2.11) 



9 



Since (2.10) holds for At by Theorem 4 of [JS], it suffices to show that an 
estimate like (2.10) holds for the function [ct 1 ^ 2 ] + \ — ct x l 2 . For simplicity 
take c = 1. We need to estimate 



\ j\[s 1/2 ] + \ ~ s 1/2 ) ds = \ J\] + \-r)r dr 

[* 1/2 1 fk l 2 f tW 1 

Y/ (k r)rdr+- ([t 1/2 ] H r)rdr. 

fcj Jk-i 2 * J[tV*\ 2 



(2.12) 



-[tV2] 



Now j k _ x {k — \—r)rdr = — |, so the first term in (2.12) is exactly 
0(t~ l l 2 ). For the second term we note that the integrand is 0(t 1 ^ 2 ) and 
the interval of integration has length at most 1, so again the contribution is 

o(rV2). □ 

The almost periodic function g is the same for An and Ad, and aside 
from the factor 1/4 it is given explicitly in [JS]. The error estimate 0(t -1 / 2 ) 
in (2.10) is somewhat worse than the estimate 0(t -3 / 4 ) for At given in [JS]. 
We may also regard (2.10) as the first term in an asymptotic expansion with 
0(t~ 1 / 2 ) replaced by sums of g k (t 1 / 2 )t~ k / 4: . For odd values of k the g k are 
almost periodic functions as given in [JS] arising from \At- The expression 
(2.12) may be written as g2(t 1 / 2 )t~ 1 / 2 + gA{t l ^ 2 )t l where gi and g± are peri- 
odic of period 1. This is easily seen because h(x) — x—[x] is such a periodic 
function, and (2.12) is a polynomial in h{t 1 / 2 ) and t x l 2 divided by t. 

Next we consider a right isosceles triangle that is half of a square of 
side length a. Then Neumann and Dirichlet eigenfunctions extend by even 
and odd reflection to eigenfunctions of the same type. Continuing the same 
notation as before with b = a, the Neumann eigenfunctions are 



c(j, k) + c(k,j) for < j < k 
and the Dirichlet eigenfunctions are 

s(j, k) — s(k,j) for < j < k. 
Lemma 2.3. For the right isosceles triangle we have 



N N (t) = -N T (t) + 



at 1 ' 2 



TV 



+ 



+ 



at 1 / 2 
V2n 



+ 



+ 



(2.13) 



(2.14) 



(2.15) 
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N D (t) 



at 1 ' 2 



+ 



at 1 ' 2 
V2n 



+ 



+ 



8 



(2.16) 



Proof. In the generic case < j < k there are eight eigenfunctions e(±j, ±k) 
and e(±k, ±j) contributing to Nt for one eigenfunction (2.13) or (2.14) con- 
tributing to Nn or Nd- The case (j,k) = (0,0) contributes the constant. 
The case k > j = contributes e(±k, 0) and e(0, ±k) to Nt, but only 



at 1 / 2 



m 



c(k, 0) + c(0, k) to A^jv and nothing to N D , leading to the term \ 

(2.15) and (2.16). When j = k > we have the eigenvalue 271 " 2 2 fc2 , with 
e(±k, ±k) contributing to Nt, only c(k, k) contributing to N D , and no con- 



tribution to N N . This leads to the term £ 



qf 1 / 2 



in (2.15) and (2.16). □ 



Thus we define the refined asymptotics 

'2 + V2 S 



N N (t) 



8tt 



t + 



Att 



at 1 ' 2 + 



(2.17) 



N D (t) 



8tt 



2W2 

47T 



at 1 / 2 + 



8' 



(2.18) 



Theorem 2.4. For the right isosceles triangle we have the estimate (2.10) 
for A N and A D 



Proof. Same as for Theorem 2.2. 



□ 



Note that 



_L + JL _L JL 
16 ~ 32 ~ 32 



as predicted for the angles |, 



7T 7T 

2' 4' 4' 



Next we consider the equilateral triangle. For simplicity we assume the 
side length is I. Six equilateral triangles tile a regular hexagon, (see figure 
2.1), and Dirichlet and Neumann eigenfunctions extended by even or odd 
reflection extend to periodic functions on the plane with respect to the lattice 
associated to the hexagonal torus T. 

We may choose a basis (\/3, 0) and (^, |) for the lattice Jif, and a basis u = 
(73' I) an< ^ = (0, |) for the dual lattice JSf*. Then the torus eigenfunctions 
are of the form 

e(k,j) = e 2n{ku+3vyx k,j E Z (2.19) 
with eigenvalue ( {k 2 + j 2 + fcj). 
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Figure 2.1 



To describe the Neumann and Dirichlet eigenfunctions on the triangle 
it is convenient to think of the dual lattice Jzf* as made up of the origin 
surrounded by concentric hexagons. In Figure 2.2 we show one such hexagon 
with a generic choice (k > j > 0) of twelve points associated to the same 
eigenvalue, together with three reflection axes of the triangle sides. 



I 




Figure 2.2 

A Neumann eigenfunction must be symmetric with respect to the three re- 
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flections, so there are two eigenfunctions associated to these lattice points, 
namely 



e(k,j) +e(-k,k + j) + e(-k-j, k) + e(—j, -k) 

+e(j, —k - j) + e(k + j, -j) 



(2.20) 



(dark points) and the same with j and k interchanged (open points). Sim- 
ilarly a Dirichlet eigenfunction must be skew-symmetric, so again we find 
two, 



e(k, j) - e(-k, k + j) + e(-k - j, fc) - e(-j, -k) 

+e(j, —k - j) - e(k + j, -j) 



(2.21) 



and the same with j and k interchanged. For the nongeneric cases we have 
(0,0) contributing to Nt and Nn but not No, (k,k) contributing six terms 
to Nt and one each to Njy and No, and (k, 0) contributing six terms to Nt 
and two terms to N N , namely e(k, 0) + e(—k, k) + e(0, —k) and e(—k, 0) + 
e(k, —k) + e(0, k) and nothing to N D . This yields 



N N (t) = -N T (t) + 



N D (t) = -N T (t) - 
o 



Thus we define the refined asymptotics 

N N (t) 

1 



-j-t 

47T 

3 

47T 



1/2 



1/2 



+ 



+ 



+ 



+ 



1 V3 3 1/2 1 
4tt 4 4tt 3' 



N D (t) = ~t-^+ 1 -, 
An A An 3 



(2.22) 
(2.23) 

(2.24) 
(2.25) 



and the analog of Theorem 2.2 is valid. Note that 1 = ^ + 5 + 5 as predicted. 



Finally, we consider the 30° — 60° — 90° triangle that is half of the equi- 
lateral triangle. For Neumann eigenfunctions that means we take Neumann 
eigenfunctions on the equilateral triangle that are symmetric with respect to 
reflection in the x-axis, and similarly for Dirichlet eigenfunctions with skew- 
symmetry. For a generic choice k > j > we get exactly one eigenfunction 
of each type by taking the sum of the two eigenfunctions of the form (2.20) 
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or (2.21), with j and k interchanged. For the nongeneric cases we have (0, 0) 
contributing to Nt and Njy but not Np, (k, k) contributing six terms to Nt, 
one term to Nn and nothing to Np, and (k, 0) contributing six terms to Nt, 
one term to N^ and nothing to Njj. Thus 



N N {t) 



—N T (t) + - 
12 w 2 



47T 




4/2 



47T 



+ 



+ 



12' 
(2.26) 



N D (t) = ±N T (t) - \ 



47T 




Thus the analog of Theorem 2.2 holds for 

N N (t) 



47T 



1 3 + yg 1/2 5 

47T 8 + Stt + 12' 



3 + V3 i/ 2 ^_ 
8tt 12' 



12' 



(2.27) 



(2.28) 



(2.29) 



Note that is the length of the perimeter of the triangle, and ^ 



Je + I + ih as P re dicted. 



3 Mixed Boundary Conditions 

In this section we look at examples of polygons with Neumann conditions 
on part of the boundary and Dirichlet conditions on the rest of the boundary. 

We begin with the simplest case: a rectangle R with Neumann conditions 
on facing edges of length a and Dirichlet conditions on the facing edges of 
length b. Just as in the pure Neumann and Dirichlet example in section 2 we 
may extend eigenfunctions to the same torus T, and so the eigenfunctions 
have the form 

3 ^ 

f(j,k) = cos7r-xsin7r-y for j > and k > 0, (3.1) 
a b 
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with eigenvalue ?Pf- + 1L ^ L - Just as in Lemma 2.1 we have 



N(t) = -N T (t) + 



1 


'at 1 / 2 ' 


1 


'bt 1 ' 2 ' 


2 


71 


" 2 


71 



1 

4 



\N T (t) + i 



it 1 ' 2 



7Y 



+ 



1 



7T 



+ 



1 

4' 



(3.2) 



and so we define the refined asymptotics 
N(t) 



ab 2a — 2b -, /9 1 

47T 47T 4 



(3.3) 



Then the analog of Theorem 2.2 holds. Note that 2a — 2b gives the differ- 
ence of the lengths of the portions of the perimeter where Neumann and 
Dirichlet boundary conditions hold, and the constant is — \ because each of 
the four vertices has mixed boundary conditions on the incident edges, with 

V»W - V(f ) = 



_j_ 

16" 



Next we consider the case where the same type of boundary conditions 
hold on one pair of opposite edges (say the ones of length b), while for the 
other pair we have Neumann on one side and Dirichlet on the other. We call 
these the NM (Neumann/mixed) and DM (Dirichlet/mixed) cases. Here we 
need a larger torus V of size 2axAb on which to extend the eigenfunctions. 
The eigenfunctions then have the form 



f(j, k) = 



COS 71 -X COS 71 

a 



sm7T-:r COS7T 

a 



y j > 0, k > (NM) 
y j > 0, k > (DM) 



(3.4) 



with eigenvalue 



ttV n 2 {2k + lf 



Ab 2 



Now we observe that N T >(t) also has eigenvalues + - f^ 1 , so N T >(t) — 
iVr(t) counts all eigenvalues of the form (3.5) with multiplicity 4 for the 
generic case j > and 2 for the case j = 0. Thus 



(3.5) 



N(t) = ±(Ni,(t)-N T (t))±± 



7Y 



t'/ 2 + 



(3.6) 
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(+ for NM and — for DM cases), since the condition ^ (2 4 ^+ 1)2 < t means 
< k < [^t 1 / 2 — |] . Thus we choose the refined asymptotics 

~ . . ab 2b , /9 
N NM (t) = irt^—t 1 ' 2 

ab 2b ^ 

and the analog of Theorem 2.2 holds. Indeed N T ,(t) - N T (t) = ™t - ^t, 
and [H 1 ' 2 + §] - ^ 1/2 = Oit- 1 ' 2 ) as in the proof of Theorem 2.2. The ex- 
planation for the coefficient of t 1//2 in (3.7) is that the sides of length a cancel 
because they have mixed boundary conditions, while the sides of length b 
add because they have like boundary conditions. There is no constant term 
because there are two vertices with mixed boundary conditions and two with 
like boundary conditions on their incident edges. 

The last case of mixed boundary conditions on R has mixed conditions 
on both pairs of opposite edges (MM). In this case we need a still larger 
torus T" of size 4a x 4b on which to extend the eigenfunctions. Then the 
eigenfunctions have the form 

j + fk+\ 



f(j, k) = costt ±j xcostt ^— ^ j y forj>0,A;>0 (3.8) 
with eigenvalue 

^ + i, V W 

4a 2 4b 2 

If we denote by T[ and T' 2 the tori of sizes 2a x 4b and 4b x 2a, then 
Nt" — N t ^ —N t ^ +Nt counts all eigenvalues of the form (3.9) with multiplicity 
4, and all eigenvalues are generic. So 

N MM (t) = \ {Nt" (t) - N n (t) - N n {t) + N T (t)) (3.10) 

Thus we define the refined asymptotics 

~ . . ab , 
N M M{t) = —t, (3.11) 

47T 

and the analog of Theorem 2.2 holds because \ (n t » (t) - N T > (t) - N T ^ (t) + N T 

^Ymm(0- Here there is no t 1//2 term in (3.11) because of cancellation of op- 
posite sides, and no constant term because of cancellation of vertex pairs. 
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Next we look at flat cylinders C obtained from the rectangles R by iden- 
tifying the opposite edges of side length b. We impose either Neumann (N), 
Dirichlet (D) or mixed (M) boundary conditions on the other pair of opposite 
edges. Then reflection produces eigenfunctions on the torus T3 of size a x 2b 
(N or D) or T' A of size ax 4b (M). The eigenfunctions are given by 



f(j, k) 



with eigenvalue 



e 2m J a x CQS7r hy 

e 2m i x sin7rfy 



J2ni L x 



a COS 7T 



4n 2 j 2 



+ 



+ 



7T 2 k 2 
b 2 

7T 2 (2fc+l) 2 



4b 2 



for j G Z, k > 
for j eZ,k>0 
for j eZ,k>0 

(N or D) 
(M). 



(N) 
(D) 
(M) 



(3.12) 



Reasoning as before we find 

N N (t) = lN n (t) + 



at 1 / 2 
2ir 



+ 



at 1 /' 2 
2n 



N D (t) = ±N T ,(t) 
{ N M (t) = l(N T ,(t)-N T ,(t)). 
We define the refined asymptotics 

N N {t) = ^t+^t 1 ' 2 



Ait 

N D (t) = ft 



N M (t) = 

and the analog of Theorem 2.2 holds. 



4-7T 

ab 
An 



A-k 

2ay-l/2 
Att L 



(3.13) 



(3.14) 



(3.15) 



Next we consider a Mobius band obtained from the rectangle R by iden- 
tifying the sides of length a with reversed orientation, with either Neumann 
or Dirichlet boundary condition on the boundary curve of length 2b. We can 
extend eigenfunctions to the torus T, and we find they have the form 



f(j, k) 



e ni i x cos ^fy 



■ j 

e m a x sin 



if k > and j + k is even (N) 
if k > and j + k is odd (D) 



(3.16) 
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with eigenvalue ^j- + 7L ^-- Consider the lattice counting function 

N e (t) = #{(j, k) : j + k is even and ^ + ^ < t} 
N (t) = #{(j, fc) : j + k is odd and ^ + ^ < i}. 

Reasoning as before we find 



(3.17) 



N N {t) = \N e {t) + 
N D (t) = lN (t) - 



at 1 / 2 



+ 



2n 

at 1 / 2 , 1 
2tt "'" 2 



(3.18) 



Now N e (t) is equal to the counting function corresponding to a torus T' h of 
area 2ab, and iV (i) = N T (t) — N T ^(t), so both \N e (t) and |iV (i) are covered 
by the results of [JS]. Thus we define the refined asymptotics 



N N {t) 
N D (t) 



ab-t _i_ 2a+l/2 

ab 
An 



t - £fV2, 



(3.19) 



with the analog of Theorem 2.2 holding. Note that the refined asymptotics of 
the Mobius band is identical to that of the cylinder, although the eigenvalue 
counting functions are not equal, and the almost-periodic function g is not 
the same. 



Next we consider the isosceles triangle with mixed boundary conditions. 
There are essentially four different cases. In the first two we consider the 
same type of boundary condition on the equal sides, and the opposite type 
on the hypotenuse. We write (ND) for Neumann on the equal sides, and (DN) 
for Dirichlet on the equal sides. In the ND case we take an odd reflection 
in the hypotenuse to obtain a pure Neumann condition on the square. If we 
compare with Lemma 2.3 where we took an even reflection, we see that the 
eigenfunctions on the square with Neumann conditions split into the mixed 
ND eigenfunctions and the pure N eigenfunctions, so 

N ND (t) + N N (t) = N { N S \t), (3.20) 

where the right side denotes the counting function for the square. Similarly 

N DN (t) + N D (t) = N { D S \t). (3.21) 
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Thus, N ND (t) is given by the difference of (2.6) (with b = a) and (2.15), so 



N ND (t) = ±N T (t) + \ 



at 1 ' 2 



7T 



+ 



at 1 ' 2 



7T 



Similarly N DN is the difference of (2.7) and (2.15), so 



N DN (t) = ±N T (t) - \ 



it 1 ' 2 



TV 



1\ 1 

+ 2 + 2 



Here the constant is V(f ) + 2 (^(|) - ^(f )) 
This leads to the choice 

^ = 8^ + H^ 



at 1 / 2 
V2n 



+ 



16 V32< 



1 

8' 



it 1 / 2 



(3.22) 



(3.23) 



(3.24) 



N DN {t) = ^t + 
and the analog of Theorem 2.4 holds. 



-2 + y/2 

An 



it 1 ' 2 



(3.25) 



In the last two cases we have mixed boundary conditions on the two equal 
sides and Neumann (MN) or Dirichlet (MD) conditions on the hypotenuse. 
By reflecting evenly or oddly in the hypotenuse we end up in the (MM) 
case for the square. As in the earlier cases the (MM) eqigenfunctions on the 
square (3.8) yield the (MN) eigenfunctions f(j,k) + f(k,j) and the (MD) 
eigenfunctions f(j, k) — f(k,j) (k ^ j). In the place of (3.20) and (3.21) we 
have 

N MN (t) + N MD (t) = N$ M (t), (3.26) 

but this does not allow us to immediately compute the two summands. How- 
ever, in the generic case j ^ k we obtain one eigenfunction in each case for 
two in the square, while in the nongeneric case j = k we obtain one in the 
(MN) case and none in the (MD) case for one in the square. Thus 



N MN (t) = ^N$ M (t) + 



V2a 1/2 
2tt 



+ 



(3.27) 
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and 



N MD (t) = ±NjS> M (t) - \ 



This leads to the choice 



V2a l/2 1 
2tt 1 + 2 



(3.28) 



N MN (t) = ^-t + 

07T 



V2 

47T 



at, 



(3.29) 



AW*) 



8tt 



47T 



at, 



(3.30) 



with the analog of Theorem 2.4 holding. Here the constant is (ip(n) — VKf )) + 

W) + Wi)-iKS)) = o. 

In our last example we consider two mixed boundary problems on the 
30° — 60° — 90° triangle with the same type of boundary condition on the 
hypotenuse and the shortest side, and the opposite type on the side that 
cuts the equilateral triangle in half that we call (ND) and (DN). In the (ND) 
case an odd reflection across the cut side produces a Neumann boundary 
condition on the equilateral triangle, so these eigenfunctions together with 
the pure Neumann eigenfunctions split up the pure Neumann eigenfunctions 
on the equilateral triangle. Thus A^jvu(t) is given by the difference of (2.22) 
and (2.26), hence 



N T (t) 



= —N T (t) 
12 w 




(3.31) 



Similarly N DN (t) is the difference of (2.23) and (2.27), hence 



N DN (t) 



—N T (t) - - 
12 TW 2 



4/2 




Att 



+ 



12' 



(3.32) 
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This leads to the choice 



N ND (t) 



1 

— —t H — 

4tt 8 8tt 



1/2 _ 



1 

12' 



(3.33) 




It does not appear that mixed boundary problems on the equilateral tri- 
angle or the other types of mixed boundary problems on the 30° — 60° — 90° 
triangle have eigenfunctions that may all be described by trigonometric poly- 
nomials. 



We consider the 2-dimensional sphere S 2 , and for simplicity we take the 
radius equal to 1. The eigenvalues of the Laplacian are given by the theory 
of spherical harmonics, with eigenvalue k(k + 1) having multiplicity 2k + 1, 
for k = 0, 1, . . . The eigenfunctions are the restriction to S 2 of homogeneous 
harmonic polynomials of degree k. We thus have 



4 The Sphere 



N{t) = k 2 for k 2 - k < t < k 2 + k, 



(4.1) 



which we simplify to 




(4.2) 



Now we define the refined asymptotics 



N(t)=t+-. 



(4.3) 



Lemma 4.1. For k 2 — k < t < k 2 + k we have 



A(t) = Y t (k 2 - (k 2 - t) 2 ) - l 



(4.4) 
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If* 1 

Proof. A(t) — - (N(s) — s) ds — -. Now for j < k we observe that 
Uo 3 
r r' ■ j rf+j 

/ (N(s) — s) ds = / (j 2 — s) ds = because j 2 — s is a linear function 

If 1 1 

on the interval varying between j and — j. Thus A(t) — - J (k 2 — s) ds — -, 
and when we evaluate the integral we obtain (4.4). □ 
Theorem 4.2. Let 

g(x)= 1 --2(x-[x+ 1 -}) 2 . (4.5) 

Then g is a continuous periodic function of period 1 with mean value zero, 
and 



A(t) = g(Jt + \) + 0(t- 1 ' 2 ). 



(4.6) 



Proof. It is clear from the definition that g is continuous and periodic, and 
if x = k + r with —\ < r < \ then g(x) = | — 2r 2 . One easily checks that 
rh i 

/ (- — 2r 2 ) dr = 0, so g has a mean value zero. □ 

If we write \Jt + \ = k + r then k 2 = t + \ - 2r^t + \ + r 2 and (A; 2 - 
t) 2 = ^ + r 2 -2r^t + iy = (\ + r 2 ) 2 -r(l + Ar 2 )^t+\ + Ar 2 (t+l 



2\2 _ r 2 



so k 2 - (k 2 - t) 2 = t(l - 4r 2 ) - yjt+\(l - 4r 2 ) + | + r 2 - (± + r 2 ) 
= i(l - 4r 2 ) - sjt + \r(l - 4r 2 ) + (4r2+3 j^ 4r2) . Thus (4.4) becomes 

A(t) = 1 - 2r 2 - r(l-4r 2 )y£ + j (4r 2 + 3)(1 - 4r 2 ) 
W 6 2 t 32t ' 1 ' ; 

Corollary 4.3. There exist continuous periodic functions g\ and g 2 such that 



t + \ 92(y/t+\ 



A{t) = g ^ f - | + ,/, ! V 'f + - | ^- ( I..8) 

More precisely, 



(x-[x + l])(4(x-[x-l]Y) 



g±(x) = and (4.9) 
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92{X) 



32 



x 



X + 



+ 3) 1-4 



x 



X + 



(4.10) 



Proof. When x 
(4.7). 



\Jt + \ we have r = x — [x + |] , so (4.8) is exactly 



□ 



Next we consider the hemisphere H, with either Neumann or Dirichlet 
boundary conditions. Reflecting an eigenfunction symmetrically or skew- 
symmetrically about the boundary produces an eigenfunction on S 2 , so the 
eigenvalues are the same k{k + 1), and the multiplicity 2k + 1 on S 2 splits 
into k + 1 for Neumann and k for Dirichlet boundary conditions on H. The 
analog of (4.1) is thus 



N N (t) 



k 2 + k 



N D (t) 



k 2 -k 



(4.11) 



For this reason we choose 



N N {t) = \t + \^t + \ + \ 
N D (t) 



\t \ J t + I + g. 



(4.12) 



We note that we could replace \Jt + \ by y/t, as that would only change 

the refined asymptotics by 0(t -1 / 2 ) and this would get absorbed into the 
remainder term in A(t). 



Theorem 4.4. Let g(t) be as in Theorem \.2. Then 



(4.13) 



for either boundary condition. 

Proof. We may write A(t) = ^Ai(t)±^A 2 (t) where Ai is the average function 
for the sphere, and 



Mt) 



If 

t Jo 



j-\js + -\ ds 
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where y s + \ = j + r with — \ < r < \. We break the integral up into the 
■ j < s < j 2 + j for j < k, and the final interval k 2 — k < s < t. 

i n r 1 * i 

j — Ws+-I(is= / 2r(j + r) <ir = -, so the contributions 

from the intervals j < k add up to ^(k — 1) = 0(t -1 / 2 ). For the final interval 
we note that the integrand is bounded and the interval has length at most 
2k, so its contribution is also 0(f) = 0(r 1/2 ). Thus A 2 (t) = 0(r 1/2 ). □ 

We note that the yJt+\ coefficient in (4.12) is exactly ^ where L = 2tc 
is the length of the boundary circle of H. The circle does not contribute to 
the constant terms because it is a geodesic in the sphere. 

Finally we consider the projective sphere PS 2 obtained by identifying 
antipodal points in S 2 . Then the eigenfunctions are exactly the spherical 
harmonics of even degree. Thus for k 2 — k < t < k 2 — k we have 

rfc-ii 

L 2 J ( k 2 +k ■£ i ■ 

„r/,\ v^/,- -,\ if k is even. , , ,\ 

fW=E«; + i)= 4, if , isodd . («4) 

Thus we define the refined asymptotics 

W) = \t + \, (4.15) 
and the analog of Theorem 4.4 holds, since the average value of the function 

{-i) k k is c^r 1 / 2 ). 



5 Lunes 

Fix an even number 2m, and consider the lune L 2m cut from the sphere by 
two great circles from pole to pole separated by the angle Let (x,y,w) 
denote the coordinates in M. 3 . We may also think of spherical harmonics as 
polynomials in (z, z, w). Eigenfunctions on the lune with Neumann boundary 
conditions extend by even reflection to the whole sphere, and so are given 
by spherical harmonics that are invariant under the rotation z h-> e~ z and 
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the reflection z (->■ z. Similarly, with Dirichlet boundary conditions and odd 
reflection the spherical harmonics muct be invariant under the rotation and 
skew-symmetric under the reflection. 



Fix a nonnegative integer N, and let F N denote the space of polynomial 
(not necessarily harmonic) that are homogeneous of degree N and satisfy the 
conditions 

f(e m z,e rn z,w) — j(z, z,w) and (5.1) 

f(z,z,w) = f(z,z,w). (5.2) 

It is clear that the functions \z\ 2 i(z m + z m ) k w e , for 2j + mk + £ = N, belong 
to Pat, and in fact give a basis for P^. Similarly if denotes the space of 
polynomials homogeneous of degree iV satisfying (5.1) and 



f(z,z,w) = -f(z,z,w) 



(5.3) 



instead of (5.2) then consists of functions of the form (z m — z m )f 1 (z,z, w) 
for /i G Pv-m (in particular Q N is {0} unless N > m). We also observe that 
A : P^ — > Pn-2 and A : — > Qn-2 and both mappings are onto. It 
follows that the multiplicities of the eigenvalue N(N + 1) in the Neumann 
and Dirichlet spectra of L 2m are given by dimP^r — dimPN-2 and dimQ^ — 
dimQ A r_ 2 respectively. It turns out to be easier to compute these differences 
of dimensions rather than dimP^r and dimQjv. 

Lemma 5.1. For any N > we have 

dim Pat - dimP^ 

while for any N > m we have 

dimQjv — dimQjv_ 2 = 



N 
m 



+ 1, 



N 
m 



(5.4) 



(5.5) 



Proof. The mapping / h-> \z\ 2 f maps Pat-2 into Pa? and is one-to-one, so 
dim Pat — dimPAr_ 2 is equal to the dimension of the complementary space, 
which is spanned by (z m + z m ) k w e for mk + £ = N. For each choice of k 
with mk < N there is a unique choice of £, namely £ = N — mk. But the 
number of choices of k with mk < N is exactly [^] + 1, proving (5.4). Then 

dirnQ^ - dimQAT-2 = dimP^_ m - dimP^-i = [^] +!=£]■ □ 



25 



Let N N (t) and N D (t) denote the eigenvalue counting functions for the 
Neumann and Dirichlet eigenfunctions on L^m- Then for k 2 — k<t<k 2 + k 
we have by Lemma 5.1 that 



fc— 1 

N N (t) = 

j=0 







m 






k- 


-i 






' f 

m 



and 



'j=m 



(5.6) 



(5.7) 



Lemma 5.2. Suppose k = p mod m and k 2 — k < t < k 2 + k. Then 

and (5.8) 



„ r . s k 2 k p(m — p) 



k 2 k p{m-p) 
ND{t) -2^-2 + ^n— 



(5.9) 



Proof. The sum in (5.6) repeats each integer from 1 up to [— ] exactly m 
times, and then [— ] +1 exactly p times. Note that k = p+m [— ] , so N N (t) = 

i ta (ta + + 1)= (ta + 1) (t [a +?)= (Hh m 

yields (5.8). Then (5.9) follows since N D (t) = N N (t) - k by (5.6) and 
(5.7). 



□ 

Now we note that the average value of the constant term as p 

varies from to m — 1 is 

Lv^ n _^ - m2 ( m - 1 ) rn(m - l)(2m - 1) 

m 2 Z-JF 1 * P ) ~ 4m 2 1W 



2m 2 



p=0 



12m 2 



1 

12 



m 



rn 



This leads to the choice of refined asymptotics 



v ; 2m 



+ 



12 



m 



+ 



6m 



(5.10) 



(+ for Nn and — for Nd)- Note that we could also replace y t + | by 
since the difference is 0(t -1 / 2 ). 
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Theorem 5.3. There exists a continuous, periodic function g of mean value 
zero such that 

A(t) = g(Jt + ±) + 0(t- 1 '*). (5.11) 
Proof. We may write D(t) = N(t) - N(t) as D^t) + D 2 (t) + D 3 (t) where 

— 2m~ {2m + 6m) ' 

D 2 (t) = ±(§-V^+P), and (5.12) 



If we form the corresponding average functions A 1 ,A 2 , A 3 it suffices to show 
(5.11) for each separately. Now the result for Ai is Theorem 4.2 multiplied 
by For A 2 it is easy to see that in fact A 2 (t) = 0(t -1 / 2 ). It remains to 
consider A 3 . 

Assume k 2 — k < t < k 2 + k for k = jm + p with 1 < p < m. We may 
write 



1 m Mnm+q) 2 + (nm+q) i rt 

4W = tEE D 3 (t)dt+- D 3 (t)dt. 

n=0 q=l Hnm+q) 2 -{nm+q) 1 J (jrn+1) 2 - (jm+1) 



It is easy to see that the last term is 0(t -1 / 2 ) because the integrand is 
bounded and the length of the interval is O^ 1 / 2 ). In each integral in the 
sum D 3 (t) is constant, so 



Mt} = iEE2(n m + 9 ,(^-l( m -I)) +0 , t -V 2) 

n=0 q=l v \ / / 



n=0 q=l 



since 



E2nm | — — ( m ) ) = 0. But the summands are 
V 2m 12 V ml J 

q=l x \ / / 

uniformly bounded and the number of terms is 0(t 1//2 ), so A 3 (t) = O^" 1 / 2 ). 

□ 
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We observe (5.10) has the predicted form. Conpared with (4.3) for the 

sphere, the leading term is reduced by the factor ^ because the area of the 

lune is times the area of the circle, and the ±| factor in the square root 

term is where 2-7T is the length of the perimeter of the lune. Finally, the 
constant is twice ^ (m — — ) , the value predicted for each of the two angels 
— of the lune. We could also obtain the analog of corollary 4.3 for the re- 
mainder term in (5.11). 

1 12 

Next we consider the half-lune L 2m obtained by slicing along the equa- 
tor. Functions satisfying Neumann or Dirichlet conditions on the equatorial 
edge extend to the entire lune by even or odd reflection; so we may study 
four different types of boundary conditions. We denote by N^, N^, 
and Np the corresponding eigenvalue counting functions, where the ± de- 
note the equatorial boundary conditions. Then N^(t) = N^(t) + N^(t) and 
No{t) = Np(t) + Np(t), and we need to understand how the spherical har- 
monics that contribute to N^(t) and Nr>(t) split into even and odd functions 
in the w variable. It is clear that this is determined by the parity of £ in the 



basis element \z\ 2j (z m + z m ) k w i . 



Lemma 5.4. (a) Assume m is even. Then 



dimP+ - dimP+_ 2 = 



dimP^ - dimP^_ 2 



dimQ+ 



dimP+ - dimP+_ 2 



dimQ+_ 2 



dim — dim ( 



"N-2 



(b) Assume m is odd. Then 



[-1 + 

LmJ 


1 if N is even 





if N is odd 





if N is even 


[-1 + 

LmJ 


1 if N is odd 


f \.m\ 


if N is even 


I o 


if N is odd 


r o 


if N is even 


l [£] 


if N is odd 




if g] is odd 



(5.13) 
(5.14) 
(5.15) 
(5.16) 



1 \in + I 

2 LmJ ~ 2 

\ i~] + ^ [~] ^ S even an d N i< s even (5.17) 

\ [m\ im\ ^ S even an d N is odd 
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if g] is odd 

dimP^r — dimP^_ 2 = < \ [^] if [^] is even and iV is even (5.18) 

i/ [^] is even and A" is odd 

[m\ ^ s even 

dimQ+ - dimQ+_ 2 = { | g] + § i/ g] is odd and JV is odd (5.19) 

^ [m] ^ S 0C ^ anC ^ ^ ^ S even 

dimQ^ - dimQ^_ 2 = { \ g] - ± i/ g] zs odd and JV zs odd (5.20) 

^ [m] ^ S °dd an< ^ N is even 

Proof. As in the proof of Lemma 5.1, we need to count the number of so- 
lutions of mk + £ = N but now with £ restricted to be even and odd. If m 
is even then £ = N — mk will have the same parity as N, and all g] + 1 
choices of k are allowed. This proves (5.13) and (5.14). A similar argument 
proves (5.15) and (5.16) since there are only [— ] choices of k. If m is odd, 
then £ will be even if A" and k have the same parity, and £ will be odd if 
A^ and k have opposite parity. If A" is even, then the even choices of k in 

< k < g] will contribute to dimP^ — dimP^_ 2 while the odd choices of 
k will contribute to dimP^ — dimP^_ 2 . If [— ] is odd there will be exactly 

1 \—] + \ of even and odd choices, while if [— 1 is even there will be \ \— 1 + 1 

2 LmJ 2 i r An 2 LmJ 

even choices and \ g] odd choices. If, on the other hand, A" is odd, it is the 
odd choices of k that contribute to dimP^ — dimP^_ 2 and the even choices of 
k that contribute to dimP^ — dimP^_ 2 . Putting this together yields (5.17) 
and (5.18). Then (5.19) and (5.20) are obtained by replacing N by N — m, 
which flips the parity of both A" and [— ] . □ 

Lemma 5.5. Suppose k = p mod m and k 2 — k < t < k 2 + k. 

(a) Assume m is even. Then 

J H + f + if pis even 

N NW - I fc 2 , /i '"! i x , , (m-p)p , m-p •, ■ , , ( 5 - 21 ) 



4m 



+ (1 + 2^+^ + ^ VP* Odd 



■ k i (m—p)p r 

+ 4 + 4^ if P is even 

- + Q-^+4? + ^ if P* odd 



N N(t) - ^ _fe2_ , /! ™ ! x , , (m-v)v . m.-n ., . , , (5.22) 

4m 

k 2 k _i_ (rn-p)p 



Nr>(t) — { k 2 /i 4m it fm-p) P P -c '• jj (5.23) 
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k I (m-p)p 

4m 4 "'" " 



4m 



4m 



i/p is even 

Assume m is odd. Then 

k 2 ( 1 1 \ . pirn — p) m — p 



(5.24) 



^(0 



4m 
4m 



4 4m 



+ 



A; 2 / 1 

+ -T + 



4m 



1 

4m 

1 



k 



k + 



4m 4m 
p{m — p) m — p 



4m 



4m 



+ h(n,p) 
— h(n,p) 



(5.25) 
(5.26) 



4 

1 

4 



, p(m — p) m — p 1 

4m y 4m 4m 2 

1 \ , p(m — p) m — p 1 , . . 

U + ^__ - + --h(n,p) 5.28 

4m/ 4m 4m 2 



where 



h{n,p) 



if n is odd 
| if n is even and p is odd 
— \ if n is even and p is even. 



(5.29) 



Proof. As in the proof of Lemma 5.2, we need to sum the expressions in 
Lemma 5.4 for TV < k — 1. Assume m is even. Write = mn + p. Then 



m) = E 



m 



+ 1 



and 



^(0 = E 



2j + l 



m 



+ 1 



The sum for has all integers from 1 to n repeated y times, and the 
integer n + 1 repeated | times if p is even and 2±I times if p is odd. Thus 



|n(n + 1) + |(n + 1) if p is even 
fn(n+l) + £±i(n + l) if pis odd. 



Substituting n = 2 we obtain (5.21). Similar reasoning shows 



Nn® 



+ 1) + |(n +1) if p is even 
m(n+ 1) + ^(n + 1) if p is odd. 
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and this reduces to (5.22). 



Similarly we have 



3< l 



2 1 

m 



and 



N D (t)= 



fc-2 



2j + l 



m 



so we just have to replace n by n— 1 in the expressions for N^(t) and N N (t). 
This yields (5.23) and (5.24). 

Next consider the case when m is odd. For N^(t) we need to sum over N 
in < N < k - 1 the values \ g] + \, \ g] + 1 or § g] according to the 
parity cases in (5.17). As before there will be m choices of iV with [^] = j 
for < j < n — 1, and p choices of N with [^] = n. So 



ra-l 



JV+(t) =m J2p + -pn 



3=0 



N 
m 



is odd} 



+ #{iV: 



AT 
m 



is even and A 7 is even} 



If n is odd then 



#{7V: 



A 7 " 
m 



is odd} = m 



n — 1 



and 

#{iV: 
In this case 



A" 
m 



is even and N is even} 



m + 1\ ( n + 1 



„. . , m(n — l)n pra m(n — 1) + (m + l)(n + 1) p 
N£(t) = 4 + y + — ^— ^ + \ 



m(n + l)n n + 1 f n + 1 



+ 



+ 
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and substituting n = we obtain 



m + 1 \ p(m — p) | m — p 



4m 



4m 



+ 



4m 



If n is even then 



TV 
m 



#{iV: 



AT 
m 



#{iV: 
is even and N is even 



is odd} = and 



i = / { m j 1 ) 1 i + P -^ if Pis odd 
I 1^)1 + 2 if Pis even. 



So 



+ , s _ m(n - l)n pn mn+(m + l)n J 2±i p odd 
^VjvW- 4 +y+ 4 +| 2 peven 

m(n + l)n n /n+l\ 1 ... . ... 

= 4 + ^ + J P + 2 ( lf P 18 odd )- 

This establishes (5.25), and (5.26) follows from N N (t) = N+(t) + N^(t) and 
(5.8). 

To compute N^(t) we note that (5.19) and (5.20) differ from (5.17) and 
(5.18) by -§. Thus N+(t) = N+(t) - (^) and N~(t) = N~(t) - (*=!), 
and this implies (5.27) and (5.28). 



□ 



Now we define the refined asymptotics 



11/ i \ i i 

4 24 V m 8 12m 



(5.30) 



— +(--—) 

4m \4 4m / 



t 



1 1 

4 4m 



11/ 1 

t H 1 m 

4 24 1 m 



1 

- - + 



1 1 

4 + 24 



m 



1 

m 



8 12m 

1 1 

+ 



11/ 1 

t+ - + — (m 

4 24 V m 



8 12m 

1 1 

+ ~ + 



8 12m 



(5.31) 
(5.32) 
(5.33) 
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Note that the coefficients of \Jt+\ may be explained because the two 

quarter circle boundary arcs have total length 7r, contributing | to and 
— I to Np : and the equatorial boundary arc has length ^, contributing ^ 

to iV^ and and — ^ to iV^ and N^. For the constant term we note that 
the top angle — always has the same type of boundary conditions on either 
side, so it contributes ^ (m — ^) , while the bottom right angles have like 
boundary conditions for and Np, so contribute | to each, and opposite 
boundary conditions for and and so contribute — | to each. 



Theorem 5.6. The analog of Theorem 5.3 holds. 

Proof. If suffices to understand the contribution of the parity dependent 
terms in Lemma 5.5. Suppose first that m is even. When we average the 
function 



<t>(j£>) = 
over 1 < p < m we obtain 



if p is even 
if p is odd 



m— p 
2m 



m — 1 

1 ^ 2j + 1 1 /m\ 2 _ 1 
m ^ 2m 2m? V ~2 ) ~8 

so the average of <p(p) — | is 0{t~ r ). 

Next suppose m is odd. This time we have to average the functions 

m-p 
MP) = ~4m~~ 

when n is odd, and 



01 (P) 



m — p \ i if p is odd 



4m \ — | if p is even 



when n is even. Now the average value of 0o is "^m"^ = \ — g^, and the 
average value of 0i is | — g^; + 4^; = § + g^; because there is exactly one more 
odd value of p than even values of p in 1 < p < m. Thus when we average 
over n we get \. □ 
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6 Cone point singularities 



The first example of a surface with cone point singularities is the flat projec- 
tive plane obtained from the unit square by identifying both pairs of opposite 
edges with reversed orientation. This yields two singular points with cone 
angle n. This example was analyzed in [JS] where it was shown that 

N{t) = \N T {t) + \± l - (6.1) 
where T is the 2x2 torus, and the choice of ± sign corresponds to the parity 



of 



t l/2 



Thus the choice 



Nit) = lt+\ (6.2) 

leads to the analog of Theorem 2.2. We note that the two cone points yield 
2 • 2^(f ) = \ for the constant in (6.2). 

Our next example is the surface of a regular tetrahedron. This is discussed 
in detail in [GKS] . Let T be the torus associated to the lattice Jzf generated 
by the vectors (2,0) and (l,\/3). This is similar to the hexagonal torus 
discussed in section 2 but it is larger, and a fundamental domain consists of 
eight equilateral triangles. There is a two-fold covering of the tetrahedron by 
T (deleting singular points), and so the eigenfunctions on T have the form 

e(k,j) = e 2ni ( ku+jv > x (6.3) 

for u = (|, ^) and v = (0,^), generators of the dual lattice j£f*, and 
(k,j) E Z 2 , with eigenvalue 

4vr 2 

— (j* + (6.4) 
while the eigenfunctions on the tetrahedron are 

cos 2n(ku + jv) ■ x = cos 2n ^7^1 + ^"^-7 + x "^j ' 

Aside from (k,j) = (0, 0), the two choices (k,j) and (—k, —j) collapse to one 
choice in (6.5), so 

N(t) = ±N T (t) + l (6.6) 
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This leads to the choice 

m = ^t + \ (6.7) 

and the analog of Theorem 2.2 holds. Note that the tetrahedron has four 
cone points with cone angles 7r, and 4 • 2*0(|) = |. 

Next we consider a half of the tetrahedron sliced so that two adjacent 
faces are cut by the perpendicular bisectors of their common edge. The 
boundary of the half-tetrahedron consists of the cut line together with a side 
edge of the face that is entirely in the surface (see Figure 6.1). 



Figure 6.1: Half-tetrahedron with boundary line in dark. 



If we impose Neumann boundary conditions then for generic choice of (k,j) 
we have eigenfunctions 



cos 2n(ku + jv) • x + cos 2ix(—ku + (k + j)v) • x, 
while with Dirichlet boundary conditions we have 

cos 2n(ku + jv) ■ x — cos 2n(—ku + (k + j)v) ■ x 



(6.8) 



(6.9) 



corresponding to (k,j) and (—k,k + j) in (6.5). However, in the nongeneric 
cases, if k — then we have one function of the form (6.8) and none of the 
form (6.9) corresponding to a single function of the form (6.5), and similarly 
for (2k, -k). Thus 



N(t) 



\N T (t) + \ 



2tt 



+ 



t l/2 

27 



+ 



(6.10) 
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for Neumann boundary conditions, and 



N(t) = \N T (t) - ± 



2tt 



+ 



jl/2 



2tt 



for Dirichlet boundary conditions. Thus we choose 



1 

4 



(6.11) 



(6.12) 



(+ for Neumann and — for Dirichlet) in order to obtain the analog of The- 
orem 2.2. Note that + 1 is the length of the boundary. To explain the 
constant term we have to remember that angles at corners have to be mea- 
sured in terms of the intrinsic geometry of the surface. Thus the boundary 
actually only has two corners, each of angle |, so these contribute 2 • ^ = |. 

Our last set of examples are the glued lunes, where we identify the bound- 
ary arcs of a lune in the same orientation. Here we will be able to allow any 
positive integer m for the lune L m , so the glued lune L m will have two cone 
points with cone angle ^. Eigenfunctions on L m are just spherical harmonics 
on the sphere satisfying rotation invariance under z i-> e^z. As in section 5, 
we may identify a basis of the space Pat of homogeneous polynomials of de- 
gree iV with this invariance as functions of the form \z\ 2 ^z mk w e or \z\ 2 ^z mk w e 
with 2j + mk + £ = N. For any fixed j and k with 2j + mk < N there is 
a unique choice of £, and there are two basis elements when k >^0 and one 
when k = 0. By considering the map / >->■ \z\ 2 f from Pat_2 to Pa? we may 
compute the difference dimP^v — dimP7v-2 by counting the basis elements in 
F N corresponding to the choice j = 0, so analogous to Lemma 5.1 we have 



dimP/v - dimP 



JV-2 



N 
m 



+ 1 



(6.13) 



Note that this is exactly the sum of (5.4) and (5.5), since the eigenspace for 
L m is exactly the sum of the Neumann and Dirichlet eigenspaces of L 2m - In 
particular, to get N(t) for L m we simply need to add (5.8) and (5.9), so if 
k = p mod m and k 2 — k < t < k 2 + k we have 



N(t) = k - + P{m - P \ 
m m 



(6.14) 
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Just as in (5.10) we choose 



N(t) = - + l(m--) + ^- (6.15) 
m 6 m 3m 

and obtain the analog of Theorem 5.3. The contribution of the two cone 
angles to the constant is 2 • 2*0(^) = |(m — ^). 



7 Sorting by symmetry types 

Suppose the surface has a finite group G of isometries. The group action 
then commutes with the Laplacian, so the eigenspaces are preserved, and we 
may sort the eigenfunctions by symmetry type, specifically the irreducible 
representations {iij} of G. This situation was discussed for a more general 
setting in [S], where it was suggested that asymptotically, the proportion of 

eigenfunctions of each symmetry type should be 777;, where dj = dim^- (it 

is well known that Yl d ] — This was proposed as a heuristic idea that 

is not universally valid (there are simple counterexamples), but should hold 
if there is a fundamental domain for the group whose boundary is relatively 
small. In the case of surfaces there is usually a fundamental domain that is 
a polygon, so its boundary is one-dimensional. Here we will be able to prove 
a more refined statement for some of the surfaces we have examined: the 
square torus, the square, the hexagonal torus and the equilateral triangle. In 
the first two cases the group is the dihedral group D4, and in the second two 
cases it is the dihedral group D 3 . 

The group -D4 has 8 elements that are generated by reflections in the diag- 
onals of the square and the reflections in the horizontal or vertical bisectors. 
There are 4 1-dimensional representations 1 ± ±, where the first ± indicates 
the symmetry type 

u(Rx) = ±u(x) (7.1) 

for diagonal reflections, and the second ± indicates the same symmetry 
equation where R is a horizontal or vertical reflection. There is also one 2- 
dimensional representation that we denote by 2. Write N±±(t)N and N 2 (t) 
for the counting function of all eigenvalues A < t that belong to each sym- 
metry type. (Note that an individual eigenspace may contain eigenfunctions 
of more than one symmetry type, and for N 2 (t) we count each eigenfunction 
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basis element separately.) 



For simplicity we choose a square of side length one, and this will repre- 
sent the square torus if we identify opposite sides, or the square itself with 
either Neumann or Dirichlet boundary conditions throughout (we nix mixed 
boundary conditions because they do not respect the D A action). So in all 
three cases we have 



N ++ (t) + N + _(t) + M + (t) + JV__(t) + N 2 (t) = N(t), 



(7.2) 



where N(t) is given by iVr(i) for the torus, and N^it) and Nu(t) in Lemma 
2.1 with a = b = 1. 

A fundamental domain for the D4 action is a right isosceles triangle of 
hypotenuse length ^ and the equal side lengths \ that makes up an eighth 
of the square. The key observation is that an eigenfunction of symmetry 
type 1 ± ±, when restricted to the triangle, satisfies Dirichlet or Neumann 
boundary conditions as illustrated in Figure 7.1, and conversely every triangle 
eigenfunction extends by the appropriate reflections to an eigenfunction on 
the torus or square. Thus N±(t) is exactly the counting function for the 
triangle given by Lemma 2.3 and (3.22-3.30) with a = \. Then we can use 
(7.2) to compute N 2 (t). 



Lemma 7.1. (a) For the square torus we have 



N ++ (t) = -N T (t) + - 



N + _(t) = -N T (t) 
N-(t) 



f l/2- 

27 
2tT 

t l/2- 

N 2 (t) 




+ 



+ 



+ 



+ 



+ 



1 

8 



+ 



-N T (t) - 



(7.3) 

(7.4) 

(7.5) 

(7.6) 
(7.7) 
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N 
1++ 



N 



D 
1+- 



N 



N 
1-+ 



9 

D 
1— 



(a) The square torus 



N 



N 



D 



N 
1++ 



N 
1+- 



1-+ 



W 

1— 



(b) The square with Neumann boundary conditions 



N 



D 

1++ 



D 

1+- 



D 

1-+ 



D 
1 — 



(c) The square with Dirichlet boundary conditions 
Figure 7.1 



(b) For the square with Neumann boundary conditions N ++ (t) is given by 
(7.3), N-+(t) is given by (7.5), 



(7.8) 



N + -(t) = \N MM {t) + ^ 


47T 


+ 2~ 


N^(t)= l -N MM {t)- l -(^ 


Vj^l/2 
47T 


1 

+ 2 


N 2 (t) = - A N T (t) - N MM (t) - 


( 


- f l/2- 

_2tt_ 


+9 



1 

4 



(7.9) 
(7.10) 



where N MM is given by (3.10) with a = \- 

(c) For the square with Dirichlet boundary conditions N + -(t) is given by 
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(7.4), N__(t) is given by (7.6), 



N+-{t)= l -N MM {t)+ l - 



N_ + {t) 



\N MM {t) - l - 




N 2 (t) = -N T (t) - N MM (t) 



Proof. Direct substitution of triangle formulas. 

We now define the refined asymptotics: 
(a) For the square torus 



++K ' 32tt I 4vr 8tt / 8 



1 

4 



(7.H) 

(7.12) 

(7.13) 

□ 



(7.14) 



N + -(t)~- 
N_ + (t) 

N-(t) = 



32vr \ 4tt 8tt / 8 



32tt I 4tt 8tt 



1 



-t 



327T 

N 2 (t) 



4tt 8tt J 8 

1 1 

8^~2 ; 



(7.15) 

(7.16) 

(7.17) 
(7.18) 



(b) For the square with Neumann boundary conditions iV ++ is given by 
(7.14), iV_ + is given by (7.16), 



N + .(t) 
N-(t) 



32vr 8tt 
32tt 8tt 



(7.19) 
(7.20) 
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N 2 (t) = — t-— t 1 ' 2 --- 
A ; 8tt 2tt 4' 



(7.21) 



(c) For the square with Dirichlet boundary conditions iV + _ (t) is given by the 

right side of (7.15), N (t) is given by the right side of (7.17), N ++ (t) is 

given by (7.19), AL+(t) is given by (7.20), 

Then the analog of Theorem 2.2 holds. 

Next we consider the case of the hexagonal torus associated to the hexagon 
in Figure 2.1. The group D 3 of order six is generated by the three reflections 
in the three bisectors. It has two 1-dimensional representations denoted 1± 
according to the symmetry (7.1) with respect to all three reflections, and a 
2-dimensional representation denoted 2. The analog of (7.2) is 

N + (t) + N_(t) + N 2 (t) = N T (t), (7.23) 

and we use this to solve for AT 2 (t). Since 1+ eigenfunctions may be identified 
with Neumann eigenfunctions on the equilateral triangle fundamental domain 
(shaded in Figure 2.1), and 1— eigenfunctions with Dirichlet eigenfunctions, 
we have N+(t) given by (2.22) and N_(t) by (2.23). This leads to 

^-hA t + T/ n + \ < 7 ' 24 > 
*-M = ST«-S tl/, + § « 7 ' 25 ' 

U) = ^t-\ (7.26) 

Finally we consider the equilateral triangle with either Neumann or Dirich- 
let boundary conditions. A fundamental domain for the action of D s is a 
30° — 60° — 90° triangle equal to a sixth of the equilateral triangle. Again, the 
eigenfunctions on the equilateral triangle with 1± symmetry correspond to 
eigenfunctions on the fundamental domain with boundary conditions shown 
in Figure 7.2. 
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N N 
1+ 1- 

(a) Neumann boundary conditions 



D 
1+ 



D 
1- 



(b) Dirichlet boundary conditions 



Figure 7.2 



Thus for Neumann boundary conditions, N T (t) is given by (2.26) and N-(t) 
is given by (3.32), while for Dirichlet boundary conditions N + (t) is given by 
(3.31) and N_(t) is given by (2.27). (Here we have to rescale because the 
triangles are smaller by a factor of 4^.) This leads to the following choices: 
(a) For Neumann boundary conditions 



~ , , 1 a/3 \/3 + 1 i /2 5 

N+(t) = ——t + t 1/2 + — 

+ V ; 4tt 24 8tt 12 

~ , , 1 a/3 -\/3 + 1 1/2 1 

iV_(t) = — — t + — 1 1 ' 2 

W 4tt 24 8tt 12 

(b) For Dirichlet boundary conditions 

~ , , 1 \/3 \/3 - 1 i 



/2 



12 



~ . , 1 \/3 \/3 + l I/O 5 

N(t) = — —t - 1 1 ' 2 H 

w 4tt 24 8tt 12 



(7.27) 

(7.28) 
(7.29) 

(7.30) 

(7.31) 
(7.32) 



Again the analog of Theorem 2.2 holds. 
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